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1. Introduction 

Kane's proposal^lj of a donor based solid state quantum computer in silicon has 
sparked a concerted effort to re-evaluate an atomistic view of impurities in doped silicon 
electronic devices. In the Kane quantum computer the phosphorus donor nuclear spins 
act as qubits, and single qubit operations are performed by applying radio frequency 
magnetic fields resonant with nuclear spin transitions. Two qubit operations are 
mediated through the electron exchange interaction. Application of voltages to metal 
gates above the spins (A-gates) and between adjacent spins ( J-gates) perturb the donor 
electron density around the nucleus, and thus the hyperfine and exchange interactions 
can be tuned with an externally applied electric field. 

We study these two gate-controlled interactions crucial for qubit operations: the 
hyperfine interaction between P nuclear spin and donor electron spin, and the exchange 
interaction between adjacent donor electrons. Here we modelled the effect of either an 
A or J-gate voltage as well as the effect of the location of the qubit in the silicon wafer 
device to determine the sensitivity of these interactions in relation to these parameters. 

In section |2l we discuss the approach we took to obtain the phosphorous donor 
ground state in the silicon wafer device. The donor wave function was expanded 
in a basis of deformed hydrogenic orbitals following Faulkner's approachj2] using an 
anisotropic effective mass Hamiltonian. To include the effect of the electric field and 
interface regions into the Hamiltonian we modeled the application of an electrostatic 
potential to the metallic gates above the qubits using TCAD,[nj and used a step potential 
to model the Si/Si02 and Si/back gate barrier. 

Section El discusses how we calculated the contact hyperfine interaction and the 
exchange interaction for the donor pair. Here we performed a Heitler-London calculation 
of the exchange coupling with the application of a J-gate potential, to study how the 
electrostatic potential enhances the exchange coupling. 

We present the numerical results for the hyperfine and exchange interaction in 
section HI and El We explore how the application of a gate voltage and the qubit 
position affects these two interactions. We can examine the selectivity of the gate 
potential by comparing the hyperfine interaction with the application of either an A or 
J-gate voltage, in E] We can also examine the connectivity between the donor pair by 
calculating the exchange splitting at varying J-gate voltage, inter-donor separation and 
donor depth. Finally we summarise our major findings in section El 

There is a great deal of attention on modelling these interactions, Kane^ makes 
a qualitative calculation using a hydrogenic approximation for the exchange coupling 
in bulk Si, without considering the electric field potential. Koiller et a/.[3E| studied 
the exchange coupling between a donor pair also in bulk Si, and the absence of an 
electric field. In their calculations they used an effective mass theory in which the 
expansion of the ground state donor electron wave function includes the Bloch states 
of the six conduction band minima. They approximated the coefficients of the Bloch 
functions using an anisotropic Kohn-Luttinger variational form for the envelope wave 
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function. Wellard et al.[9l have extended these calculations to remove some of their 
approximations. They obtained the donor electron wave function and bare exchange 
coupling at zero J-gate bias, in order to study the fast exchange oscillations with respect 
to fabrication strategies. 

Fang et a/. [IT] calculated the donor electron wave function using the spherical 
effective mass approximation. They modeled the J-gate potential qualitatively as a 1-D 
parabolic well with its minimum located in the middle of the two donor sites, but did not 
consider the boundary effects of the silicon host geometry in their calculation. In their 
work they used an unrestricted Hartree-Fock method with a generalised valence bond 
wave function to study the two-electron system and calculate the exchange coupling. 
Parisoli et al. have calculated the effect of the J-gate potential, interface regions and 
donor separation using a spherical effective mass Hamiltonian. We extended this work 
to include the anisotropy of the effective masses in Si into the Hamiltonian. 

We study the effect of application of a 3-D electrostatic potential to the metallic 
gates above the qubits, and boundary effects of the silicon oxide layer and back gate on 
the donor electron wave function. We included the anisotropy of the effective masses, the 
P impurity potential, electric field and interface potentials into the Hamiltonian. We 
calculated the contact hyperfine interaction and exchange coupling for varying qubit 
separation, qubit depth and gate voltage. We aim to provide relevant information for 
experimental engineering of these devices and highlight the significance of environmental 
factors other than the gate potential which may perturb the donor electron wave 
function. 

2. Faulkner's Method with the applied electric field and silicon host 



Using the method outlined previously, |1] we use an anisotropic effective mass 
Hamiltonian, Ho, for the donor in bulk Si and zero field: 



where e = 11.4 is the dielectric constant, and 7 = m±/m\\ = 0.2079. Here we are using 
atomic units, where the unit of length, as = h^e/m^_e^ = 31. 7A and the unit of energy, 
Eb = m^e^/2h^e^ = 19.94meV. 

Following Faulkner's approach we expanded the donor electron wave function in 
a basis of 91 deformed hydrogenic orbitals: 



where ipnim{x,y, z,a) = Rni{a,r)Yim{0, (f)), are the normalised hydrogenic orbitals, 
Cnim are the expansion coefficients for our basis, a is the effective Bohr radius in 
the X, y directions, and f3 is an adjustable parameter which gives the effective Bohr 



potential 




(1) 
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Figure 1. Schematic design parameters implemented in TCAD to model the Kane 
computer architecture. 



radius b in the z direction. Equation ([Q) was solved variationally using a basis of 91 
deformed hydrogenic orbitals for the donor in zero field, to give a ground state energy 
E = -31.23meV, and effective Bohr radii: a = 23.8lA and b = 13.68A.|1] 

To accommodate the effect of the applied field and the boundaries on the donor 
electron wave function it is necessary to use more than one simple bulk ground state 
wave function to describe the envelope function. The method we used is advantageous 
because we expand the envelope wave function in a basis of deformed hydrogenic orbitals 
which have the flexibility to distort with the applied fields. We use the zero field effective 
Bohr radii and diagonalise the single donor electron Hamiltonian including the electric 
field. 

To include the effect of an electric field and the silicon host, we constructed an 
additional Hamiltonian matrix. Hi, with its elements given by: 

{n'l'm'\Hi\nlm) 

xi)ni^{x,y,^z,a), (3) 

where Kiecd/, z) is the electric field potential generated from TCAD, and here we also 
add an additional term to model the Si02 layer and the back gate as a step function with 



'- / dx^ip*^,i,^,{x,y,J-z,a)Veieciy,z) 
1 J V 7 
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height 3.25eV. jH El 112] The lateral edges of the silicon lattice were assumed to extend 
infinitely in the ^/-direction, and the potential in 2-D from TCAD is assumed to have a 
"thickness" in the third dimension (x) of 1/im. Figure ^ shows the 2-D device scheme 
implemented in TCAD used to model the application of voltages to the A or J-gate 
above qubit, Qi, the metallic gates were modelled as thin wires in the x-direction. 

The new Hamiltonian H = Hq + Hi was then diagonalised to find the perturbed 
single donor electron ground state for each particular gate voltage and qubit position. 
The location of the interfaces in figure Q splits the degeneracy of the two conduction 
band minima along the 2;-axis relative to the other four along the x and y-axis, in the 
lower Ai ground state in zero electric field. [Z] We expect that at the shallow donor 
depths we consider, the greatest restriction on the donor electron will be the interface 
regions. Here we expand around the conduction band minimum in the 2;-direction from 
Qi to the silicon oxide layer. This was done so that the smaller effective Bohr radius, 
6, would be in the direction toward the silicon oxide and back gate barrier. Using this 
convention the donor wave function is lower in energy since there is less overlap of the 
wave function into the interface regions. 



3. Calculation of the hyperfine interaction coupling and exchange splitting 

For the Si:P quantum computer to be feasible, quantum operations have to be able to be 
applied selectively to particular nuclear spins, and connectivity between nuclear spins 
via electron-mediated coupling must be established. To achieve both these goals it is 
necessary to study the degree of selectivity and connectivity that can be controlled by 
applying electric fields to metal gates above (A-gates) and adjacent (J-gates) to spins. 
Furthermore, it is shown in this paper that the qubit location in the device in relation 
to each other (inter donor separation) and to the gates (donor depth below the silicon 
oxide barrier), also has a significant infiuence on the donor electron wave function. 



3.1. Calculation of the contact hyperfine interaction 

Since we use effective mass theory, instead of calculating the contact hyperfine coupling, 
^4(1^), directly we calculate the relative shift in AiV) with the potential applied and 
assume this shift will be similar to those of the true wave function. jHl E] Thus we need 
to calculate: 

where A{Q)/h = 28.76MHz is determined for ^^P in silicon from experimental data. [TBI IT] 
and \E'(V,r) are the donor envelope wave functions calculated by our method. 

The contact hyperfine interaction was calculated for the varying J-gate voltage, 
inter donor separation R, and donor depth below the silicon oxide layer, and compared 
with our previous results for similar calculations at varying A-gate voltages. |3] 
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3.2. Calculation of the exchange splitting for an impurity pair 

In this section we employ a Heitler-London (H-L) treatment of the two electron donor 
pair wave function, using the two single donor ground state wave functions perturbed 
by the electric field as our basis. Since the donor ions are generally well separated in the 
silicon wafer device we can justify using H-L theory to describe the two electron system 
as the symmetrised and anti-symmetrised products of the single donor orbitals at each 
qubit (\I''5^(r) and "^^^{r)) calculated with the electric field applied. The singlet and 
triplet impurity donor pair wave functions are given by:j 



T 



1 



where 

^orbit 

J2{1±S\ 

^Qi(j.)^*Q2(r-R)c/rl 



^Qi(j.i)^^2(r2 - R) ± ^^Hr2)^^'(ri - R) 



(5) 



Here ^'^^(r, V) and ^^^(^r, V) are the single wave functions calculated using our basis of 
deformed hydrogenic orbitals, and diagonalising the Hamiltonian for the varying voltages 
at the J-gate and qubit position. We observe that '^'^^{x,y,z) = vl>Q2(^ 

the 

donor wave functions on adjacent nuclei are mirror images about the y-axis when a 
voltage is applied to the J-gate (see figure [T} . 

To calculate the exchange splitting between the ground singlet and triplet states 
for an impurity pair of donors in silicon we use the H-L formula: J3] 

J{R) =Et-Es 

= {^T\H2e\^T) - {^s\H2e\'^s) 

= {s'Ko - K,) , (6) 

2 2 2 2 



where: 

H2e = - V^„j^(ri) - V^„j^(r2) 



kil |r2| Iri — R| |r2 — R| 



+ ] 7 + Velec{ri) + V;iec(r2), 



Ki = J ^*«i(r2)^*'5^(ri - R)^^i(ri)^^2(r2 - R)Qdrldr. 



e 



2 2 2 

|r2-ri| |ri-R| |r2| 



4. Results obtained varying gate voltage and inter donor separation 

To demonstrate the effect of J-gate voltage and inter donor separation on the donor 
electron ground state we calculated the perturbed single electron donor ground states 
as a function of these external factors. The calculations in this section were obtained 
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Figure 2. Contact hyperfine interaction at varying gate voltage and inter donor 
separation, for a donor depth of 20nm. 

at a donor depth of 20nm. Once the perturbed ground state under the apphed field 
was obtained we calculated the contact hyperfine interaction and the exchange splitting 
(using (jD) and (0) respectively), for the impurity donor pair, to optimise and determine 
the experimental conditions needed to control the nuclear spins coupling to the donor 
electron spin, via the hyperfine interaction, and to other nuclei via the electron-mediated 
exchange interaction. 

So far we have only considered the effect of the A or J-gate independently. The 
smaller inter donor distances [R < 14nm) are only possible if the gate dimensions can 
be reduced to prevent overlapping gates. For this work we have only considered the 
impact of the J-gate without reference to the A-gate. In this initial study we aim to 
give insight into, and identify the relevant factors that contribute to the hyperfine and 
exchange coupling, which need to be studied more in depth. 

Results for the contact hyperfine interaction 

Figure El shows the magnitude of the contact hyperfine interaction calculated for varying 
qubit separation and gate voltage. As the donor electron density decreases at the P 
nucleus so does the contact hyperfine interaction. These results refiect the trend that 
as the qubit moves away from the J-gate, the donor electron wave function has more 
freedom to move towards the J-gate and distort greater. 

Figure El shows an example of the donor ground state wave functions of Qi and Q2 
for an applied voltage of 1.0 V at the J-gate, for two inter donor separations, R = 14 and 
20nm. From the relative magnitudes of the ground state electron densities of the two 
qubits we can see that the ground state wave functions for R = 20nm have perturbed 
more towards the J-gate voltage. This figure demonstrates how the electron density 
is perturbed greater for larger inter donor separations, (which also implies a greater 
distance of the qubits from the J-gate) at large positive gate voltages. 
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Figure 3. Ground state electron densities of Qi and Q2, in y-direction at z = 0, a; = 
for an inter donor separation of i? = 14 and 20nm, a donor depth of d — 20nm and a 
voltage of 1.0 V at the J-gate. 



We observe in figure El for certain negative gate voltages that the contact hyperfine 
interaction A{V) ~ 0, which indicates that the donor wave function has distorted 
completely away from the nucleus. For R < 14nm, and V < —0.5V, the electron is no 
longer bound to the nucleus, and disperses completely away from the applied voltage. 
Similarly for 16 < i? < 20nm, and V < —0.6V, the electron is no longer bound to the 
nucleus. This is also reflected in figure lU where we see an abrupt change in the donor 
wave functions going from i? = 14 to 16nm for a voltage of -0.5 V at the J-gate. For 
R = 14nm the electron density at the nucleus is close to zero and the wave function 
disperses from the negative voltage in all directions. In contrast the wave function for 
R = 16nm is still bound to the nucleus and only perturbed slightly by the applied 
negative voltage. 

The effect of the gate voltage on the donor electron depends on the distance of the 
qubit from the gate.jU The donor depth of 20nm and inter donor separations 

considered in this section {R < 20nm) means that the qubits are situated at relatively 
short distances from the J-gate. So for positive voltages the electron transfer to the 
gate with increasing J-gate voltage is gradual. (1^ However for negative J-gate voltages, 
there is an abrupt change in the electron density at the nucleus for critical negative 
voltages where the electron is no longer bound to the nucleus. We find depending on 
the distance from the gate and the magnitude of the negative gate potential, the electron 
transfer is either gradual or abrupt. 

4-2. Results for the exchange splitting 

In order to compare our results with previous work[71 |H1 El UHl HH IE]; we evaluated the 
zero field exchange interaction, J{R), in bulk Si, for varying inter donor separation. Our 
results using effective Bohr radii, a = 2.381nm and b = 1.368nm, are in close agreement 
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Figure 4. Ground state electron densities of Qi and Q2, in y-direction at z = 0, a; = 
for an inter donor separation of i? = 14 and 16nm, a donor depth oi d — 20nm and a 
voltage of -0.5 V at the J-gate. Note the different vertical scales of l^'p in (a) and (b). 



to the calculations of Fang et a/.^I] The zero field exchange splitting calculated here 
is higher than other reported theoretical values using H-L theory, [71 |H1 IHl UDI because 
we chose the larger Bohr radius, a = 2.381nm, to be along the inter donor axis, and 
hence the exchange splitting is larger using this convention. The larger Bohr radius was 
chosen to be along the inter donor axis so that it would also be towards the positive 
J-gate potential, and the smaller Bohr radius in the direction towards the interfaces. 

Until calculations are performed which include the effects of the interfaces on the 
donor wave function, it is hard to verify whether a higher exchange energy would 
in fact be expected because of the decreased probability of penetration of the donor 
wave function into the interface regions. Koiller et a/.0 calculated the exchange 
coupling in uniaxially strained Si in the presence of interfaces, and also found that these 
environmental influences could affect the exchange coupling significantly. They found 
that the F±{z) envelopes were favoured energetically, because the smaller effective Bohr 
radius in the z-direction guarantees less signiflcant penetration of the wave function into 
the barrier regions. In this paper we are trying to model the effects of the electric field 
potential and Si host geometry on the donor wave function, to investigate the variation 
of the exchange splitting with the applied voltage, rather than the absolute values of 
J{R). 

Figure El presents our results for the exchange coupling as a function of inter donor 
separation and positive J-gate voltage. We observe that the exchange coupling increases 
as the J-gate voltage increases as expected, since the applied field draws the electrons 
closer together. At a voltage of l.OV the donor electron wave function is perturbed the 
greatest, and the exchange coupling is significantly higher at this voltage for every inter 
donor separation. 

At all voltages lower than l.OV, the exchange coupling decreases as R increases 
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Figure 5. Calculated exchange coupling as a function of inter donor distance and 
positive J-gate voltage, (a) is the results for V < 0.7V and (b) for V = l.OV, for a 
donor depth of 20nm. 
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Figure 6. Calculated exchange coupling as a function of inter donor distance and 
negative J-gate voltage, for a donor depth of 20nni. 



as expected, but in figure Efb) for a voltage of l.OV and for R = 20nm, the exchange 
couphng actually increases slightly. This is because at large inter donor separations the 
donor is further from the J-gate, and thus is more attracted to the potential well at the 
J-gate as the gate voltage increases sufficiently. Hence the overlap between the adjacent 
donor electron orbitals is slightly greater, even if the inter donor separation is higher. 

One of the advantages of the Kane quantum computer is the ability to turn on 
or off the coupling between the different qubits. Figure IHl presents our results for the 
exchange coupling as a function of inter donor separation and negative J-gate voltage. 
The exchange coupling decreases as the negative applied potential decreases. When the 
applied negative voltage is large enough, the electron is greatly distorted away from 
the nucleus. We have seen this in figure IH^a) that at a J-gate voltage of -0.5V and 
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R = 14nm, the donor wave functions for Qi and Q2 have perturbed away from the 
apphed voltage in opposite directions. In this case we have effectively turned off the 
coupling between the adjacent qubits, as the overlap between the two electron densities 
is almost zero. 

Fang et al. ^Tj used a spherical effective mass Hamiltonian and modeled the J-gate 
potential using a one dimensional parabolic well, and in their work they do not consider 
the effect of interfaces. In their calculations of the exchange splitting they considered 
relatively large inter donor separations (> 16nm) and a smaller electric field potential at 
the J-gate relative to the TCAD cross-sectional potential we obtain between the donors. 
This made it hard to compare to our work, but as a rough estimate we can compare our 
results for a voltage of 0.2 V at the J-gate, where the TCAD potential at the mid-point 
between the two donors, approximately 0.02 V, is nearly equal to the potential at the 
minimum well used by Fang et a/.|ll^ for /i = 0.6. 

We observe that our calculations are up to an order of magnitude lower than 
those calculated by Fang et al. This discrepancy is probably because the electric 
field potentials used by our method and Fang et al. are not equivalent, as we obtain 
the electric field potential within the whole device, whereas they use a simplified 1-D 
potential. We also compared our results with Parisoli et al. ^01 IH] who used a spherical 
effective mass Hamiltonian, with effective mass, m^, ^ 0.29mo and effective Bohr radius, 
a ~ 2nm. The results for the exchange coupling agreed qualitatively (to within an order 
of magnitude) and predicted the same trend in variation of the exchange coupling with 
voltage. Again our results were consistently higher because of the larger effective Bohr 
radius along the inter donor axis. 

The application of a voltage to the J-gate, and the magnitude of the qubit 
separation, can be used to control the strength of the exchange coupling of the donor 
pair. At lower voltages the most significant factor influencing the exchange coupling is 
the inter donor separation. As the positive voltage is increased, the donor electron wave 
function is perturbed greater as the donor is moved further away from the gate, and 
the exchange coupling can be enhanced even with large qubit separations. At smaller 
inter donor separations the donor electron is more affected by negative voltages. For 
voltages lower than a certain critical value, the donor electron is completely transferred 
away from the nucleus and the negative voltage, and the exchange coupling decreases 
almost to zero. 

5. Results obtained varying the donor depth 

We observed the effect that the donor depth, d, below the silicon oxide layer has on 
the donor ground state perturbed by a gate voltage. Here we ran calculations for 
5 < d < 45nm. 
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Figure 7. Contact hyperfine interaction at varying donor depths and J-gate voltage, 
with inter donor separation R — 20nm. 



5.1. Results for the contact hyperfine interaction 

Figure [7| show our results for the contact hyperfine interaction at varying donor depth 
and gate voltage. Here we observe similar trends in the variation of the hyperfine 
interaction as in the previous section. We see that as d increases, and hence the 
distance from the J-gate also increases, we see a cross-over behaviour where the donor 
wave function is perturbed greater for larger donor depths at positive voltages above a 
critical value. Also for d > 25nm we see in figure [Z(b) there is an abrupt decrease in 
the electron density, defining an ionisation voltage at these donor depths. This process 
of ionisation has been reported previously. |31 

For d > 25nm and large enough positive voltages, the electron has perturbed or 
ionised completely to the gate. Figures 01 and |H1 show the contrast in the donor wave 
function for different donor depths. In figure |H1 at c? = 30nm, the electron is perturbed 
almost completely away from the nucleus towards the applied voltage. While in figure El 
at d = 20nm, the donor wave function is only slightly perturbed from the zero field 
ground state by the applied voltage. 

The basis we are using for the donor electron wave function consists only of bound 
states, which is a good approximation for the smaller gate voltages, as the electron is 
still bound to the nucleus. To model the ionisation process at larger gate voltages more 
accurately, a more rigorous approach would be to include the delocalised conduction 
band states in the basis as well. 

Figure Ufa) and (b) shows the comparison between the electron ground state 
probability density in the j/z-plane, obtained for a voltage of l.OV applied to the A and 
J-gate respectively, for a donor depth of 30nm and inter donor separation of R =20nm. 
For an A-gate voltage the donor wave function is symmetric in y and only perturbs 
toward the A-gate in the z direction. In comparison, when a J-gate voltage is applied, 
the wave function can distort in both the y and 2;-directions. Unfortunately as the donor 
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Figure 8. Ground state electron density in yz-plane for donor depth at 30nm with 
l.OV at the A-gate in (a) and J-gate with R = 20nm for (b). In both plots Qi is located 
at the origin, and we have included the y = symmetry line in the contour plot, to 
highlight the difference in the electron density for an applied A or J-gate voltage. 

depth becomes greater, selectivity may be lost, and a voltage applied at either the A or 
J-gate will cause the same change in the contact hyperfine interaction. 

5.2. Results for the exchange splitting 

Figure El shows the variation of the exchange coupling with donor depth and voltage 
for two inter donor separations, R = 14 and 20nm. It is evident that the depth of the 
donor influences the degree to which the electron is perturbed by the gate voltage, and 
hence will also affect the strength of the exchange coupling. 

For small d the electron is only slightly perturbed by the positive gate voltage, as 
the P nucleus and the gate voltage are so "strongly coupled" [IHI that the effect of the 
gate voltage is only to further stabilise the donor electron. So we see for d = 5nm that 
although the exchange coupling has increased significantly from the zero field coupling, 
it is still not as strong as the coupling for d = 10 and 20nm. 

The exchange coupling for d = 10 and 20nm are similar for equal inter donor 
separations, and the effect of the magnitude of d is not so pronounced. Here we can 
see that for a donor depth of lOnm the exchange coupling is enhanced the most by 
the applied voltage. One of the reasons for this may be that for d = lOnm the donor 
wave functions predominantly move toward the apphed voltage at the J-gate in the y- 
direction along the inter donor axis and thus the exchange coupling is enhanced further, 
whereas for d = 20nm the wave functions can perturb in both the y and ^-directions 
toward the J-gate. 
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Figure 9. Calculated exchange coupling as a function of donor depth and J-gate 
voltage, with i?=14nm in (a), and i?=20nm in (b). 



6. Conclusions and prospects for achieving silicon-based quantum 
computation 

In this work we have studied the P donor wave function perturbed by an electric field 
and the Si host geometry, and the two interactions fundamental to the Kane quantum 
computer: the hyperfine and exchange interactions. We have studied the effect of 
varying several experimental parameters: the gate voltage, inter donor separation, and 
donor depth in order to fine tune the hyperfine and exchange interactions. 

The results presented highlight the significance of not only the gate potential in 
affecting the donor electron wave function, but also the position of the qubits in the 
device. One of the critical discoveries was that the inter donor separation is not the only 
relevant factor in determining the strength of the exchange coupling, the proximity of 
the qubit to the gate is also important in determining the degree to which the electron 
exchange interaction can be enhanced by the applied voltage. 

In the absence of an electric field, only the inter donor separation is instrumental 
in determining the strength of the exchange coupling, and as R increases the exchange 
coupling decreases. However, when a large positive voltage is applied at the J-gate, 
either a gradual transference of the donor electron density occurs for dopants close to 
the gate, and the exchange coupling is enhanced proportionally. Or if the electron is 
ionised by the gate voltage the exchange coupling can be enhanced considerably even 
for quite large inter donor separations and donor depths. So both of these competing 
influences must be considered in modeling the strength of the exchange coupling. 

For the parameters we studied (ic. R < 20nm) and negative J-gate voltages, 
V < —0.6V, the electrons at Qi and Q2 disperse away from the applied negative voltage 
at the J-gate in opposite directions. In this case the overlap between the two electron 
densities is almost zero, and the exchange coupling also decreases almost to zero, and 
we have effectively turned off the coupling between adjacent qubits. 
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Future developments in our laboratory is concentrating on confirming these results 
using a more rigorous evaluation of the exchange coupling. However this initial work 
provides valuable insight into the environmental influences which play a vital role 
in determining the sensitivity of the donor electron wave function to the applied 
electric field and hence the sensitivity also of the exchange coupling and the hyperfine 
interaction. 
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